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Abstract 

Spin transport driven by the temperature gradient in ferromagnetic metals is studied based on a 
microscopic theory. It is shown that the temperature gradient works as an effective field equivalent 
to the electric field as for both the spin current generation and the spin relaxation torque. The 
thermally driven contribution of the spin current and the relaxation torque are thus proportional 
to VT and V 2 T, respectively. 
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I. INTRODUCTION 



Thermoelectric effects such as the Seebeck and the Nernst effects have been studied for 
more than a hundred years, and are applied to various devices like thermocouples. The 
effects have been successfully explained by phenomenological and microscopic theories^, as 
for the charge transports. 

Recently, spintronics, which aims at the control of the electron spins, is attracting special 
attention from the viewpoints of the fundamental science and application. Of particular 
importance in the spintronics is the spin current. The spin current is generated by applying 
the electric field^ or by use of the magnetization dynamics via the spin pumping effect^. 
Detection of the spin current can be carried out electrically by use of the inverse spin Hall 
effect^, which converts the spin current into the charge current or electric voltage using the 
spn-orbit interaction. In 2008, another method to creat the spin current was discovered by 
Uchida et al., namely, the spin Seebeck effect^. They have applied a temperature gradient 
to a ferromagnetic metal (permalloy) under the condition that no charge current flows. By 
attaching a Pt contact on the permalloy in the perpendicular direction, they found that there 
appears a finite voltage across the Pt lead. Since no charge current flows in the permalloy, 
Uchida et al. concluded that the voltage is the result of the inverse spin Hall effect due to 
the spin current that is induced by the temperature gradient in the permalloy. Uchida et 
al. thus demonstrated that the spin current can be induced by the temperature gradient 
similarly to the Seebeck effect for the charge. 

Thermal effects on magnetic domain walls such as the eddy current induced at the domain 
wall due to the Nernst-Ettingshausen effect were discussed by Berger-~— . Very recently, 
theoretical studies on the thermally-driven spintronics phenomena have been intensively 
carried out^~— . Thermal spin-transfer torque was discussed by Hatami et at^, and its 
inverse effect was argued by Kovalev et al.— . An unified description of magnetic, electric, 
thermal and mechanical forces was presented by Bauer et al.— . 

In the conventional (charge) Seebeck effect, the effective electric field E is induced pro- 
portional to the temperature gradient, VT, as E = SVT, where S is called the Seebeck 
coefficient. As general argument indicates that the Seebeck coefficient of free electrons at 
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low temperatures is written by a energy derivative of the electric conductivity 0B(e) as^ 

2 da B (e) 

S = —{k B ) 2 T * 
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(1) 



where \x is the chemical potential. In ferromagnets, the exisitence of the charge current 
indicates also that of the spin current, since the conduction electrons are spin polarized. 
Defining the spin current in a uniform ferromagnet as j s = j + — j_ , where j± denotes the 
current carried by the electron with spin ±. Treating the two spin channels as independent, 
thermally induced spin current reads from Eq. ([TJ 

^2 



VT, (2) 



where <tb,± represents the conductivit for the spin ± electron. We define the spin Seebeck 
coefficient S s as j s = cr s S s VT, where a s = o"b,+ — cb,- is the spin cunductivity. It then reads 



3e yB) E±(±)^b,±(6) 



(3) 
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The aim of this paper is to derive this expression on a microscoic model, and to extend 
the argument to a general case with inhomogeneous magnetization, and to study the spin 
relaxation torque. 

Most crucial feature of the spin current when compared to the charge current is the 
violation of conservation law in solids. The spin density s and the spin current density j s 
thereby satisfy the continuity equation 

s a + V-j: = T a , (4) 

where a — x, y, z is the spin index and T is the spin relaxation torque resulting in the non- 
conservation of the spin. In metals, the dominant origin of T is the spin-orbit interaction. 
The relaxation torque has been treated by introducing a phenomenological spin chemical 
potential and the spin relaxation time^^. The relaxation torque plays essential roles in 
spintronics phenomena such as the current-induced magnetization switching^ and the in- 
vserse spin Hall effects^. Thus microscopic study of the torque is urgent and important. 
The spin relaxation torque induced electrically was recently studied microscopically^, but 
the thermal contribution has not been explored. 
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The aim of this paper is to theoretically study the spin transport induced by the tem- 
perature gradient. The temperature gradient is modeled by considering a system made up 
of subsystems having different temperatures and chemical potentials. Each subsystem is 
assumed to be in local equilibrium. The electron transport is studied by introducing the 
electron hopping between subsystems. For spin current, we also take into account the in- 
homogenuity of the magnetization, up to the first order in the spatial derivative. The spin 
relaxation torque is studied in the homogeneous magnetization case and in the presence of 
the spin-orbit interaction due to random impurities. In the context of charge transport, the 
temperature gradient has been known to act as an effective electric fieldi. We will show the 
equivalence of VT and the electric field holds also in the spin transports and in the spin 
relaxation phenomena. 

Although the spin-orbit interaction is essential in studying the spin relaxation torque, we 
will neglect it in studying the spin current, since we are interested in how the temperature 
gradient acts as an driving force on the dominant spin current but not in deriving the full 
transport equation. In fact, the spin-orbit correction to the spin current has been known 
to have the same dependence on the driving field as the contribution without the spin-orbit 
interaction^. 

II. MODEL 

We model the temperature gradient by considering a discretized model consisting of 
the systems labeled by n = 1,2- ■•. Each system n is assumed to be in local thermal 
equilibrium at temperature T n and chemical potential fi n (Fig. [1]). (In the end, we will 
take the continuum limit, assumeing that the temperature gradient is not very large. ) 
Without losing generality, we assume that systems are placed on a cubic lattice with equal 
distance d. The conduction electrons in each subsystem are represented by plane waves 
whose wave vectors k are approximated to take any value. The magnetization direction of 
each system, n n , is assumed to be uniform within the system but is different for different n. 
The Hamiltonian of the systems when isolated is given as 




(5) 



n 
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where M is the spin spitting energy due to the magnetization. The electron operator is 
represented by a two-component field, c n u = (c n fc+,c n fc_), where ± represents the spin. To 
describe the magnetization n n dependent on n, a gauge transform in the spin space that 
diagonalize the exchange interaction is useful. This transform is carried out as 



(6) 



where a n k is a new electron operator in the gauge-transformed frame and U n is a 2 x 2 
unitary matrix, given as 



U n = m n -cr 



[sin — cos (p n , sin — sin <p n , cos — 



with 



being the polar coordinates representing n 
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(7) 
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FIG. 1: The discrete model we consider, made up of subsystem labeled by l,m,n- ■ ■ connected 
by leads. Each subsystem n is assumed to be in the local equilibrium at temperature T n and the 
chemical potential /i n , and have a uniform magentization n n . The center coordinate of the system 
n is represented by X n , and the spacing of the systems is d (i.e., \X n — X m \ = d for a neighbouring 
pair. The electron hopping occurs on the lead (shown by obals) , between sites and Rnrn • 

By the gauge transformation, the Hamiltonian of the subsystems when isolated becomes 

Hq = £k± a ^ nk± a nk±, (8) 

n fe± 
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where €k± = J^- =F -M", and ± is the spin index. 

The subsystems are connected by leads, where the electron hopping occurs. The coor- 
dinate in the lead in a system n where the hopping to a neighbouring system m occurs is 
represented by R nm . The hopping Hamiltonian reads (in the real space representation) 

Ht ^ ^ ^ ^ ^ (^m±(-^ini)^i±(-^nm) C^n± ( Rnm)Cm± ( Rmn) ) i (9) 

(nm) RnmRmn ± 

where (nm) denotes a pair of neighbouring systems. After the gauge transform, it reads 

Ht ^ ^ ^ ^ t(o) n (R nm )U nm Cl m (Il mn ) -\- m ( Rmn) U mn (X n ( R nm ) ) , (10) 

(nm) RmnRnm 

where £/ nm = U ] n U m . 

In terms of a n electron, the spin density of the n system is written as s n = (c^crcn) = 
(alUlcrU n a n y We represent the charge and spin currents through the junction as and If 
(i = x,y,z and a = x,y,z are the spatial and spin direction, respectively). In the present 
junction model, the spin (charge) current at system n is calculated by estimating the time 
derivative of the spin (charge) density, which reads 

^re(-^wi) ^ ^ ^ ^ ^ ( < ^n(-^nm)^n ( ^" U n U nm (l m (R mri ) (l* rn (R mri )U mn U^ l O' U n Q, n (R nm )'^ . 
'in Ji mrl 

(11) 

The current flowing in direction % is thus given by 

1% ( n ) = 2 ^ ^ ^ tr[C/^(T ^([/^(^^(JZu,,,, ^^-(^^(iinm, ■Rmn)^mn)]x m =X„±iii) 

ran Rnm- 

(12) 

where i is a unit vector along z-axis and m is the label of the system neighbouring the 
system n (i.e., X m = X n ± id), and < denotes the lesser component. Estimating the 
Green's functions to the lowest (second) order in t, we obtain 

w = -^E E (±M£^n 

i fimil Rnm R' m n R'nm 



x {U nm g m {Rmni R mn )U mn g n (R nmj R nm ) g n (R nm , R nm )U nrn g m (R mn , -K mn )f^mn)]x TO =x m ±id' 

(13) 

where g m is the free Green's function of the system m on the Keldysh contour. Without 
losing generality, we choose n n as along z direction, i.e., U n = 1. Since we are eventually 
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interested in the continuum limit, we consider the case where the difference between n m 
and n n is small. The rotation matrix U nm is expressed by a spin gauge field A nm defined 
as U nm = e lAnm ' cr = C7^ n . Explicitely, U nm = m n ■ m m + ier ■ (m n x m m ) and thus 
A nm = (m n x m m ). The current then reads 

7 » = \ E(±) [y- J2 MM"" {9nk.e lA --°F mk ,(g^ - gUJe^— 
± J 71 kk' 

_l_ P l'n a - n r U^tim'ff „a -iA„ m -<T iA„ m -tr r -iA„ m '(r p / a r \ 

T" r nk\y n kuj ynkujjt ilmk'u^ e iJmk'uj^ r nk\y nkuJ iJnku}) 

— e fmk'{ymk>LU ~ 9mk'uj) e 9 n kuj f \X m =X n ±id 

l^ fc | 2 trK(F mfc , - F nfc )Im[^ fe JIm[^ fc ,J 

kk' 

± ^ k^ 

+ F nfc Im(^)[cr /3 , ^ fc , w ]-K, ^J^Jm^J-^, F mk ,lm(g^ kl J]g^ kui }] Xm=Xn±id 

= i& + «?, (14) 

where 

i fc , ;fc = t gifc'-R^g-ifc-R™ ( 15 ) 

and /q and 5/ a represent the contribution without the gauge field and the linear order 
contribution, respectively. We neglect the higher order contribution in the gauge field, since 
we consider a slowly varying spin texture. The geometry of the lead is reflected in the 
amplitude t k ',k- The Fermi distribution function is represented by matrix 

, /m(efc+) 

F mk = ( | , (16) 

f m (e k - ' 



fmi^ka) = g 7- — y (17) 



where (/3 m = (k B T m ) x ) 

Retarded Green's function g r mkuJ is a 2 x 2 matrix in spin space with each component defined 

as 

9 ™ k ™ = fa-e^ + ir)' (18) 
where rj represents an infinitesimal positive (or the inverse lifetime if disordered). 



III. UNIFORM MAGNETIZATION 



Let us first consider the contribution Iff, the current when the magnetization of the whole 
system is uniform. Explicitely writing the spin index, the currents read 



± J 71 kk' a 

/dcu 
^ ^ |tfc'fc| 2 ^o-(/m(efc'a) - /n(e fc a))Im[^ fc(T JIm[^ fc 



k'crtui X m —JC n ±id 



' cruj\X m =X n ±idi 



± " kk' 

(19) 

and Iq = Iq = 0. We consider the case of an infinitesimal difference of the temperature and 
the chemical potential of the two adjacent systems, n and m. Defining T rn = T + A, m T and 
A*m = A* + A m /i (T and fi are the temperature and the chemical potential for the electron of 
the system n), we expand the physical quantities up to the linear order in the differences. 
The difference of the Fermi distribution functions for n and m is written as 

(/m(efev) - fn(e ka )) = (/ n (e fe ' CT ) - f n (e ka )) + f'(e k ' a ) ( A m fi + (e fe v - n)—^- J , (20) 

where 

/V.) = ~{ — -. (21) 

4 cosh |(e fe / CT - /i) 

Let us first consider the conventional Seebeck effect, i.e., the charge part. The charge 
current reads 

_ 2 ^ frfw^i, ,2 1 \ \ 

6V 2 J 2tt £ 1 fc ' fcl 4- ( W - e fc .) 2 + ^ ( W - 6 fcV ) 2 + r/ 2 cosh 2 f (e fe , CT - 
E( ± )f A ^+( e ^-/ i )%^) • ( 22 ) 



X 

-j- V / X m =X n ±id 

The w-integration is carried out as 

2 



(itu 1 1 77 1 

rf I = — I h c.c. 

27r (uj - e ka ) 2 + rf (u - e k > a ) 2 + rf 2 (e feCT - e fe / CT ) V (tka - tk'o + 2ir)) 

(23) 



(efc CT - e fc / CT ) 2 + 4r? 2 
and thus 

J ° 0j = ^2 S |tfc ' fc|2 ? (^"^) 2 + V cosh 2 f(e fc , ff - /x) 

xE(±)(A m/ i+(e fcV -/i)%^ N ) • (24) 
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We consider now the continuum limit by taking d — > 0. This is allowed when d C A, where 
A is the spatial length scale the temperature and the chemical potential varies significantly 
, i.e., A ~ O ({-jr) 1 , (tt ) J ■ The summation over m is carried out easily as (choosing 
X n = 0) 



~ ^2{±)A ml J > \ Xm=Xn±id = ^H+di) - /i(0) - (/x(-dt) - /i(0))] = dVifi. (25) 
We thus see that 

j° = -d{GVn/e + G T VT) , (26) 

where the conductance G and the thermal conductance Gt as functions of the chemical 
potential are given as 

G(u) = 2e 2 VV H± 

w * (efcCT " ^ + 4 ^ cosh2 Ifa" ~ A 

M (^-^) 2 + ^cosh 2 f (efc , CT -^)' (27) 

respectively. 

The conductance is written by use of the Boltzmann conductivity <tb as G = a-^A/d 
where A is the area of each system. Let us here switch to the current density, defined as 
j = I/A. Then Eq. (126]) reduces to 

j = -(o- B V/Ve + <TrVT). (28) 

The conductivity is given in terms of the spin resolved conductivity <7b,o- (cr = ± is the spin 
index) as 

dea B ^e)—ff -. (29) 

-aM cosh |(e-/i) 

where 

'■-^^Si^W 1 '-" 1 (30) 



At low temperatures, cosh 2^ e ^ = ~ /-0 an d we reproduce 0"b = J2a a B,a- Defining 
x = |(e — /i), the thermal conductivity is written as 

ar = k ^EQ^ dx ^ aB -° u+lf - (31) 
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At low temperature, /3/j, ^> 1, we can expand the integrand with respect to x to obtain the 
well-known relation 



hn'T 7T 2 



(IT 



e 3 



E 



da 



B,cr 



de 



<Jb x 



k R 2 T7l 2 ^ ^ 



€=)! 



e=fi 



e 3 T.a^M ' 

From the above result, the spin current density, jg = Iq/A, is easily obtained as j'q 
(j s V/i/e + <t Sj tVT, where 

/3/4 



(32) 



<7, 



E 

<T=± 



dea B)<r (e)- 
ct m cosh |(e-/i) 



(33) 



and 



0"s,T = &s X 



(34) 



e 3 <7<?B,<t(J*) 

We therefore see that the temperature gradient acts on the spin current as an effective 
electric field in agreement with naive guess. When the magnetization is along n, the above 
result of spin current becomes 



j Q = -n> s V/i/e + a s , T VT). 



(35) 



The Seebeck coefficient is defined as the ratio 

cr de 



s 



e=fi 



°b e 3 Y^a a ^,M 
We define the spin Seebeck coefficient as 



(36) 



a SjT k B 2 T ir 2 



5* 



de 



e=/i 



(J, 



e 3 E ff ffff B,,(/i) 

Therefore our model reproduces the relations (JTj) (J2j) obtained by classical argument^. 



(37) 



IV. EXAMPLES 

The explicite expressions of a-B.a , S and S s depend on the detail of the hopping on the 
lead. We here present results for three typical cases. 
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A. Point- like lead 



We first consider a case of point-like lead. The coefficient tk<k then becomes a constant, 
tk'k — t (since we can choose R nm — Rmn — 0), and the summations over k' and k become 
independent. In this case, 

a B)tr (e) = 2irt 2 e 2 j "ov^fc ~ e Wa ) = 27rt 2 e 2 j(v ) 2 (e + aM), (38) 
k' 

where u = ls the three-dimensional density of states divided by y/e (e is the enery 

and V = Ad). (We note that in taking the pointlike limit, t 2 d(vo) 2 /A needs to be kept a 
constant, since the combination t 2 d(i> Q ) 2 /A gives the physical conductance as seen in the 
above equation.) We then obtain 

a T 7T 2 k B 2 T 



ob 3 efi 

while 



(39) 



*s,T = 0, (40) 

since d(JB ^^^ is a constant independent of the spin. Thus spin Seebeck coefficient vanishes if 
the lead is point-like and if the conduction electron's energy is the free electron type, k 2 . 

B. Two-dimensional interfaces 

If the junction in the discretized model is a thin plane with electron scattering, the 
electron hopping conserves the wave vector perpendicular to the junction (which we denote 
kj_) but not the component along the junction (we choose the junction along x axis). The 
wave vectors before and after the hopping are thus written as k = (k, k±) and k' = (k', k±), 
respectively, where k and k' are independent. When the continuum limit is taken in this 
interface model, the summation over the wave vectors in Eq. ( 1301) is carried out to obtain 



o- B , CT (e) = 2nt 2 e 2 j(vi ) 1) ) 2 vi ) 2) V7T^M, (41) 

where v^p = -0=^ and = ^ are the coefficients in the one- and two-dimensional density 
of states, respectively. Therefore 

a. T tt 2 k B 2 T a T 



(T, 



6 ev //i 2 - M 2 <tb' 



(42) 
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C. Free electron limit 



Our model can also describe the free electron limit. To describe the free electron with 
mass m, the parameter t and t^y of the tight-binding Hamiltonian need to be replaced 
by t ^ ^2 (using tcos(k x d) = || + const.) and t kk > = tk x d5 kk > (i.e., t\ k , = -J^S^), 
respectively. (We choose the transport as along x direction.) We then reproduce from Eq. 
fl30|) the Boltzmann conducitivity 



e 2 n rT [e)T 



*V(<0 = - J ^> ( 43 ) 
mr 

where n a (e) = ^j^^ (k(e + aM) = ^/2m(e + aMj). The spin Seebeck coefficient then 
reads (from Eq. ([57]) ) 

s = k B 2 T n 2 ( /i + M)l-(/i-M)^ 
e 2 (^ + M)§ - Gu-M)f ' 

V. GAUGE FIELD (SPIN TEXTURE) CONTRIBUTION 

We consider here the contribution linear in the gauge field, 5 If, in the adiabatic limit, 
namely, the spin variation is assumed to be small so that the momentum transfer due to 
the gauge field is neglected compared with the conduction electron's momentum. The linear 
contribution in Eq. (114j) reads 

5If{n) = -e£(±) /^El^| 2 ^ m tr[an^Ja^^Im(^ fc ,J] 

± J kk> 

± ^ fcfc' 

+ a Q [F nfe Im(^)[a /3 , <4 fe ,J - K, <4fcd^nJm(^ fc J]] Xm=Xn± ^ (45) 
By use of [a Q , ^ fc J = -2z E 7 e Q7 ^ 7 E CT °"^fcu,<r> the s P in P art (« = z, V, z ) reads 

«?(n) = + (6j? + - (46) 
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where A is the junction area and the coefficients are 
2 c f duj 

<H,P = 5^(±) h^Yl l^'kl 2 A nm a 'Unka ~ fmk'a' )Im(^JIm(^ fcw ) X „ 

± kk' acr' 

^iJ3 = T / 7T~ ^ Kfe ; fcl 2 ^nm ^ q " q "7»fc f 7lm(^ fcci;CT ^ fc , CJCT ,)x m =X n ±id 

± ^ fcfc' crcr' 

± fcfc' crcr' 

(47) 

The coefficient M 6 ' 5 -' represents the equilibrium spin current and represents the driven 
contribution. 

By writing f mk , a , = f k , a , + 5f k > a > (f k > a > = f nk i a > and 8f k > a > = 
f'i^k'a 1 ) (A m /i + (e k > a i — the summation over the spatial directions is carried out 

by expanding the gauge field and the chemical potential as (we choose X n = 0) 



\ J2(±)Ai m A m »\ Xm=Xn±id = \ £(±)(m(0) x m(±di)f(fi(±di) - /x(0)) 

d 3 (V 2 if i)Af + o(d 4 ,VA), (48) 



rj / j \ 1 Tim I -^-m— ^nl'« o 

/ ± ^ ± 



where A? = ~ (m x Vim)' 3 is the spin gauge field in the continuum limits. We have neglected 
the contribution containing the derivative of A?, since it corresponds to the second order 
contribution with respect to the gauge field, which we neglect in 51. The equilibrium and 
the linear contributions are given as (we here suppress the index n in the Green's functions) 

Oi t p = 2A^{a E {Vi) 2 fi/e + a T {Vi) 2 T) 

b^ = 2A^b E iy l Y^/e + b T iy l YT) 

b<$ = 2A*j<*>, (49) 



13 



where 



and 



a E 
ax 

b E 
b*(e) 

»(eq) _ 



3 e 



a ^7 



d 3 , 

" 2e T E - 



»?|*fc'fc| 



fcfc'"-' 



E a6 - 

3 e 



Edb a 
a ^7 



- 2e T E - 



2 



/ l^fe'fel (Cfccr — Cfccr' 



r^(e - e ka ), 



2e ~A E ^^'fel^M^L^feW'^ 



(50) 



(51) 



represents the equilibrium current. 

The result ( 1461) is for the magnetization at n is along the 2 axis (i.e., U n = 1). A general 
case with magnetization along n is obtained by applying an unitary transformation defined 
by a matrix V a p = 2m a m l3 — 5 a @. By use of the identities^ 



{2m a m p - S a(3 )(S^ - 5p z 5 lz )Aj 



(2m a m (S - 5 a p)e^ z A] = --V*n 



-(n x Vin) a 
1. 



(52) 



the final result of the gauge field contribution to the spin current density (5jf = 5 If /A) is 
given by 



5jf = -[(n x ViTi) a j^ + (V i n) a (a B V^/e + a T V 2 T) 
+ (n x V i n) a (^V 2 /i/e + 6 T V 2 T). 



(53) 



To summarize the results of the spin current, the total charge current in the system is 
given by Eq. ff28l) as 



Ji = Jo.i 



ia B V ifi/e + arViT), 



(54) 
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and the spin current reads (Eqs. (l35l)(!53l) ) 



3s,i — 3o,i + 

= n a {a s V if i/e + a^V^T) - (Vin) a (a E V 2 fi/e + a T V 2 T) 
- (n x V,n) a (j^ + ^V 2 /i/e + 6 T V 2 T). 



(55) 



The spin current driven by the field and the spin texture (spin gauge field) was calculated 
here by assuming adiabatic condition. This is justified when the electron mean free path 
t satisfies I 3> d, namely, either in the strongly disordered or in the weak hopping regime 
where I oc 1/t 2 is large. 

VI. SPIN RELAXATION TORQUE 

In this section, we calculate the spin relaxation torque induced by the temperature gra- 
dient by including the spin-orbit interaction. We consider an uniform magnetization case 
and neglect the spin gauge field. The leads connecting subsystems are assumed here to be 
point-like for simplicity, namely, = t. The spin-orbit interaction is 



The spin-orbit potential v so is assumed to arise from random impurities and to depend on 
the spin direction (k). The impurity scattering is treated in the standard manner—. 

By deriving the continuity equation for the spin density, the dominant spin relaxation 
torque in z direction acting in the system n is found to be^ 




ijk 



(56) 




(57) 



ijkl 



It is calculated including the hopping to other subsystems at the second order as 



T z (n) = --m so A so 2 t 2 ^^ / —k 2 {yye zM tx[a^g nyl ^g nku g mkl(jiJ g nh 



fh/ m kk'k" 




< 



(58) 



Taking the lesser component, we obtain 



m hh'h" J rr 




X lm (9nk",-a,J(fn(^) - fm(e k 'a) JM^iifcaJM^ow ) , 



(59) 
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FIG. 2: The Feynman diagram describing the dominant contribution to the spin relaxation torque 
induced by the temperature gradient. v so represents the spin-orbit interaction, and t denotes the 
electron hopping to other subsystems m with different temperature T m . 

where the distribution function of subsystem n is given by (j3 m = (/c^X^) -1 ) 

fm( e ka) = g 7—— r—r - (60) 

The summation over the wave vectors are calculated by use of contour integrals, and we 
obtain 

32 r°° i 

T z {n) = — — ix 2 m 2 n &0 \ S0 2 t 2 ijl V] V] gt„ / dx — (A/i m + 2k B AT m x) 

" ™ J-Q(u-M) cosh X 



m a 



fl + (7 M+^) (fi-aM+^-Y (61) 



where A/x m = fi m — fi, AT m = T m — T (T and \i are the temperature and the chemical 
potential of the system n), x = f (e^/ — crM — /i). Considering low temperatures, i.e., 
/?(// — M) ^> 1, we can expand the integrand with respect to x and obtain 



T» = - Q X 2 m 2 n so \ s 2 t 2 ul £ £ ar CT (/i 2 - M 2 ) 2 fa - aM)" 1 



1/2 
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m a 



f , 7T 2 , „ . _ 7u - crM \ , , 

x^.+ T * B TAT. 5 £-_ y j (62) 

This is the result of a discrete model. We now take the continuum limit by replacing 
Afx m with X m ■ V/i + • Xm > i ^ m -3 VjVj/i, and similar expression for AT m (we have chosen 
X n = 0). Considering the rotationally symmetric system with equal separation d between 
the small local equilibrium subsystems, X m = ±dr, where r represents three unit vectors in 
the three spatial directions, we obtain A/i m = <i 2 V 2 /i. The relaxation torque is therefore 
obtained as 

T z = 7sV 2 /i/e + 7 tV 2 T, (63) 
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where 

IE = -^Wn so \jevld> (/i 2 - M 2 ) 2 <rr. (A* - oM )~ 1/2 (64) 

cr 

7t = -p7r 4 ^ 2 Tm 2 n so A so 2 t 2 z/ 3 c/ 2 (/i 2 - M 2 ) ^ ar CT (/i - aM)~ 1/2 (7/i - aM) (65) 

cr 

The relaxation torque arises thus from the second order derivatives of /i and T . The result 
for fi here confirmes the result of Ref.— in a discretized model. (Unlike Ref.—, Eq. ( 163]) is 
symmetric with respect to the spatial direction, but this would be an artifact of the present 
model, which assumes that the electron hopping occurs on pointlike leads.) We see that, as 
is expected, the temperature gradient VT is equivalent to the electric field in the context of 
the spin relaxation. 

Equation fl65l) indicates that when an uniform temperature gradient is applied to a fer- 
romagnet, the spin relaxation torque is zero. Therefore, the spin current driven by homoge- 
neous temperature gradient in the spin Seebeck system^ is spatially uniform without decay, 
and would be consistent with the experimental observation of the inverse spin Hall signal 
over the sample of millimeter size. We stress here that a term proportional to the spatial 
coordinate introduced without ground in Ref.— does not exist in the transport equation. 

For understanding the experimental result of the thermally induced inverse spin Hall 
effect, the present analysis needs to be extended to incorporate the spin-charge conversion 
due to the spin-orbit interaction, which will be carried out in the forthcoming paper. 



VII. CONCLUSION 



To conclude, we have studied the spin current and the spin relaxation torque driven by 
the temperature gradient microscopically by considering a continuum limit of a discretized 
model. We have shown that the temperature gradient acts as the effective electric field and 
drives spin current. In the uniform magnetization case, the spin Seebeck coefficient is given 
by 

£±(±) d 



(66) 



e 3 E±(±Vb,±(/u 

where ctb,± is the Boltzamnn conductivity for the spin ± electron and \i is the chemical 
potential. 

17 



When the magnetization n is nonuniform, spin current components polarized along n x 
Vn and Vn are induced by the temperature gradient (Eq. fl55l) ). We have also calculated 
the spin relaxation torque and found that it is proportional to V 2 T. Since the relaxation 
torque induced by the electric field has been shown to be proportional to V • E, we see 
that the temperature gradient VT acts as the effective electric field in the context of the 
relaxation torque, too. 

We have thus demonstrated the equivalence of the temperature gradient and the electric 
field in the spin transport. We note, however, that quantitatively these two fields lead to 
different results since the ratio of the coefficients such as <t t /<j b and Jt/ie ar e not always 
equal. 
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